function Phi = features(state, RBFs)
    Phi = zeros(length(RBFs), 1);
    for i = 1:length(RBFs)
       Phi(i) = computeRBF(state, RBFs(i).centre, RBFs(i).width);
    end
end